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Abstract 

The distribution of the initial very short-time displacements of 
a single particle is considered for a class of classical systems with 
Gaussian initial velocity distributions and arbitrary initial particle 
positions. A very brief sketch is given of a rather intricate and lengthy 
proof that for this class of systems the n-th order cumulants behave 
as t 2n for all n > 2, rather than as t n . We also briefly discuss some 
physical consequences for liquids. 

1 Introduction 

This paper is meant to give a pedagogical sketch of the proof of a theorem 
for the distribution of initial very short time single particle displacements in 
terms of cumulants for a class of classical systems with smooth potentials, 
which can be in or out of equilibrium. This theorem is based on a time 
expansion around t = 0. 

The theorem has a number of physical applications [T], including incoher- 
ent neutron scattering in equilibrium systems E], heterogeneous dynamics 
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in supercooled glass-formers 4 and the recently developed Green's function 
approach to transport on picosecond time scales, when the system is far from 
equilibrium jSJ |H| . 

In order to formulate the theorems we first introduce the moments and 
cumulants of the distribution of single particle displacements. 

2 Moments and Cumulants 

In a classical system of N particles with positions and velocities Vj, we 
can consider the displacement Arj(t) = r;(t) — 1^(0) of an individual particle 
i in a time t. For simplicity, we will only consider here the displacement 
Axi(t) = Ari(t) • x of particle 1 in the x direction. 

Given the time t and an initial point in phase space, Axi(t) is unique. 
However, we will consider an ensemble of initial conditions, so that each 
initial point has a certain probability associated with it. In fact, the theorem 
discussed here is applicable to systems whose initial phase space distribution 
is of the form 



Here and represent the collection of the coordinates = (ri, T2, ■ ■ ■ , rjv) 
and the velocities = (vi, v 2 , . . . , v^v) of all N particles in the system, re- 
spectively, Ui = (vj)t=o for each particle i — 1,2, . . . , N and f(r N ) is a well 
behaved, i.e., normalized, but otherwise arbitrary probability distribution 
function of the r N at i = 0. 1 We will take Uj = in the following for 

In eq. Q a multi-component mixture would have different rrn for particles i, depending 
on to which component they belong, while /3, allows the particle to belong to a component 
with a different temperature at t = than other particles may have. 
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simplicity. 

Because of the ensemble distribution of initial phase space points, Ax\ (t) 
is not fixed but has a probability distribution. To characterize this probability 
distribution function, one may determine its moments 

H n {t) = ([A Xl (t)r) (2) 

The first moment (n = 1) is the average of Ax\(t), the second moment 
is the average of its square, etc. In may be noted that for small times t, 
Axi(t) ~ Vi x t, whence one expects 

fin(t) CC t H . (3) 

The moments may alternatively be calculated from /i n (t) = d n F s (k, t) / d(ik) n \k = Q 
using the moment generating function 

F s (k,t) = $> n (0^f- (4) 
= ( e ikAx ^). (5) 

Here we used eq. (J2J in the last equality. The quantity on the right hand 
side of eq. (0) is precisely the self-part of the intermediate scattering function 
(the Fourier transform of the Van Hove self-correlation function) and can be 
measured through incoherent neutron scattering |2|. Note that we need in 
principle all /i n (t) for n from zero to infinity to find F s (k,t). 

Sometimes however, it is known that the first two moments suffice, namely 
for an ideal gas and for a perfectly harmonic system, in which cases F s 
and the Van Hove self-correlation function are Gaussian. This may still be 
approximately true in less ideal situations, and in fact, for short times (such 
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as the ones that we are interested in here) one expects the system to behave 
almost like an ideal gas. So it makes sense to try to expand around a Gaussian 
ClS db zeroth order approximation. 

In such near-Gaussian cases, the so-called cumulants K n are a more con- 
venient set of parameters to work with than the moments|7j. The cumulants 
can be found from the cumulant generating function log F s (k,t), i.e., 

Mt) = ^ logjFUM) — (6) 

so that 



fc=0 



F s (k,t) = exp 



Ln=l 



(7) 



Here we will write for convenience the n-th cumulant K n (t) = ((Ax^)), 
and stress that [n] is not a power but an index indicating the order of the 
derivative of the generating function; for n — 1 the superscript will be omit- 
ted. Equating the right hand sides of eqs. (J1J) and (J7j), one can derive the 
first few cumulants in terms of the moments: 

Ki(t) = 

«a(*) = Vh(t) — 

K 3 (t) = f i 3 (t)-3 f i2(t)fi 1 (t) + 2i4(t), etc. (8) 

In general K n (t) is composed of products of moments of degree m, by 
dividing its index n in all possible ways into a sum of integers m < n, so that 
their sum m equals n. From eq. (jSJ), one would then expect: 

K n {t) cx t n . (9) 
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3 The Theorem 



SchofieldjH] and Sears found by straightforward but complicated calcula- 
tions for the initial short time behavior of the K n (t) the following results in 
equilibrium (where all odd indexed cumulants vanish): 



instead of the expected behavior K n (t) = 0{t n ) suggested by eq. ©• 

Our theorem is a generalization of Schofield and Sears' suggestive results 
for the first three non- vanishing cumulants for neutron scatterings in equilib- 
rium to the general class of systems, characterized by the initial distribution 
function ((TJ). 

The theorem states that if a) the interparticle and external forces on the 
particles are smooth and independent of their velocities, and b) the velocities 
are Gaussian distributed and independent of the initial at t = 0, then 



where the c n depend on the forces but not on t. 

This result implies hidden correlations in the moments fj, n (t), which 
are not eliminated by the simple moment expansion (jSJ) of the cumulants 

K n (t). 

4 Sketch of Proof of Theorem 

The proof of the theorem proceeds as follows, where we give only some 
of the most important steps. The full proof is presented in ref. 1. 



« 2 (t) = 0{t 2 )- K A (t) = 0(t 8 ); « 6 (t) = 0(t 12 ) 



(10) 




(11) 
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1. We consider systems of iV particles % = 1, . . . , N in d = 3 with Af = 
3N degrees of freedom. It is convenient for the formulation of the proof 
to associate with each degree of freedom of the particles i, (generalized) 
positions and velocities Vi respectively, but now with i = 1, . . . ,Af. Thus 
T\ x -> n,r ly -> r 2 ,r lz -> r 3 ,v lx = v u . . . , transforming r^.v^ to r N ,v N . 
The equations of motion are then for % — 1, . . . , M: 



Ti = Vi 



v^F^^/m^a^.t) (12) 



2. Expand Ari(i) in powers of t around t — 



3. Here the coefficients of t m are polynomials in t>i: 
P 2 , ^ = ai (^, t ) = 0(1 ); 

P 4 - - 5 ^ = ... + gg s - s - w = 0((^)') etc. (14) 
^. Therefore, the time expansion of Ari(t) can be written in the form: 

An(t) = f^Pj(v M )t j = P l {v M )t + P 2 (v u )t 2 + . .. + P n (v M )t n + 0(t n+1 ) (15) 
i=i 

where for j > 1, Pj is a polynomial of degree j — 2 in . 

5. Then the expansion of K n (t) up to t 2 "^ 1 can be written in the form: 

^(t) = «Ar> ] (t)» 




77-1 ! • • • n. 



(16) 



^ 1 =2 «3 0'-2)<ni 



with a correction of 0(t 2n ). In this equation, the semicolons on the right 
hand side of eq. (fTHj) indicate that the Pj^ inside the ((. . .)) are not to be 
multiplied, since they are elements of a cumulant. 



eq. (fTB*|) vanish, so that only the 0(t 2n ) remains and K n (t) = 0(t 2n ) (for 
n > 2). The proof of this result is based crucially on the Gaussian properties 
of the velocities, so that (vf n ) = (2n — l)\\(vf) n obtains. 

The theorem can be generalized to the displacement of a single particle in 
different directions in a d- dimensional space, as well as to the displacements 
of different particles, and can also be applied to multi-component mixtures^. 



The theorem has a number of physical applications, mostly pertaining to 
liquids, e.g. 1) It provides a well ordered short time expansion of the Van 
Hove self-correlation function relevant for incoherent neutron scattering in 
equilibrium systems 2J; 2) The cumulants are connected to the non-Gaussian 
parameters a n introduced by Rahman and Nijboer|3] and which are used as 
indicators of dynamical heterogeneity in supercooled glass-formers |3]; and 3) 
It provides a well ordered short time expansion of the Green's functions for far 
from equilibrium (mass, momentum and energy) transport on the picosecond 
time scale 1310]. For a more extensive account of these applications we refer 
to ref. 1. 



6. One can prove then that all n terms from t n to t 



2n-l 



in the sum in 



5 Physical Applications 
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